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Making use of the Stochastic Vacuum Model and the gluon gauge-invariant two-point correla- 
tion function, determined by numerical simulation on the lattice in both quenched approximation 
and full QCD, we calculate the elementary (quark-quark) scattering amplitudes in the momentum 
transfer space and at asymptotic energies. Our main conclusions are the following: (1) the am- 
plitudes decrease monotonically as the momentum transfer increases; (2) the decreasing is faster 
when going from quenched approximation to full QCD (with decreasing quark masses) and this 
effect is associated with the increase of the correlation lengths; (3) dynamical fermions generate two 
components in the amplitude at small momentum transfer and the transition between them occurs 
at momentum transfer near 1 GeV^. We also obtain analytical parametrizations for the elementary 
amplitudes, that are suitable for phenomenological uses, and discuss the effects of extrapolations 
from the physical regions investigated in the lattice. 
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I. INTRODUCTION 



Soft hadronic scattering, characterized by long distance phenomena, is one of the great challenges in high-energy 
physics. The difficulty arises from the fact that perturbative QCD can not be applied to these processes and, presently, 
we do not know even how to calculate elastic hadron-hadron scattering amplitudes from a pure non-perturbative QCD 
formalism. However, progresses have been achieved through the approach introduced by LandshofF and Nachtmann 
0, developed by Nachtmann and connected with the Stochastic Vacuum Model (SVM) 3]. In that formalism 
the low frequencies contributions in the functional integral of QCD are described in terms of a stochastic process, 
by means of a cluster expansion. The most general form of the lowest cluster is the gauge invariant two-point field 
strength correlator 0, |^ 
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-[5^(zp(5^cr - ZaStyp) + dy{z„5pp - Zp5^„)][l - n)Di{z'^ /a )}, (1) 

where z — x — y \s the two-point distance, a is a characteristic correlation length, k a constant, < FF > the gluon 
condensate and Nc the number of colours (C,D = 1, ...,N^ — 1). The two scalar functions D and Di describe the 
correlations and they play a central role in the application of the SVM to high energy scattering Q ■ Once one has 
information about D and Di, the SVM leads to the determination of the elementary quark-quark scattering amplitude, 
which constitutes important input for models aimed to construct hadronic amplitudes. Numerical determinations of 
the above correlation functions, in limited interval of physical distances, exist from lattice QCD in both quenched 
approximation (absence of fermions) 0, 0] and full QCD (dynamical fermions included) ■ 

In a previous paper, we determined the elementary amplitudes from lattice QCD in the quenched approximation 
0, using as framework the SVM. In this communication, we apply the same procedure, now taking into account the 
lattice results in full QCD. Our main goal is to investigate the differences in the elementary amplitudes associated with 
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quenched theory and full QCD and also the effect of different bare quark masses. In addition, we obtain analytical 
parametrizations for the amplitudes that are suitable for phcnomcnological uses, and discuss in some detail the region 
of validity of all the results. 

The manuscript is organized as follows. In Sec. II we recall the main formulas related with the elementary amplitudes 
in the context of the SVM and in Sec. Ill we review the parametrizations for the correlators from numerical simulations 
on the lattice. In Sec. IV we present the results for the elementary amplitudes from full QCD with different quark 
masses and discuss the similarities and differences with our previous result in the quenched approximation. The 
conclusions and some final remarks are the contents of Sec. V. 



II. STOCHASTIC VACUUM MODEL 



In this Section we briefly review the main steps of the calculation scheme that allows the determination of the 
elementary amplitudes through the SVM 0, H. We refer the reader to for more details concerning specific 
calculation. 

The elementary amplitude / in the momentum transfer space may be expressed in terms of the elementary profile 
7 in the impact parameter space, through the symmetrical two-dimensional Fourier transform 

POO 

f{q^) = / bdbJoiqbhib), (2) 
Jo 

where is the momentum transfer, b the impact parameter and Jq is a Bessel function. 

In the Nachtmann approach Q , the study of the elementary scattering is based on the amplitude of quarks moving 
on lightlike paths in an external field, picking up an eikonal phase in traveling through the nonperturbative QCD 
vacuum. In order to have gauge invariant Dirac's wave function solutions a Wilson loop is proposed to represent each 
quark. In this context the no-colour exchange parton-parton (loop-loop) amplitude can be written as 

J^{lr[Pe •> loapl —l^j.^'Pg J !oop2 ~ IJ): 

where () means the functional integration over the gluon fields (the integrations are over the respective loop areas), 
and u; is a common reference point from which the integrations are performed. 

This expression is simplified in the Kramer and Dosch description by taking the Wilson loops on the light-cone. In 
the SVM the leading order contribution to the amplitude is given by Q 

i{b) = m'{b), (3) 

where t] is a constant depending on normalizations and 



<^{b) = 9 J j da^,i,d(TpaTr{F^^{x]w)Fp„{y]w)). 

Here {g'^F^i,{x; w)Fp^{y; w)) is the Minkowski version of the gluon correlator. 

After a two-dimensional integration, e{b) may be expressed in terms of the correlation functions in (1) by Q 



e{b)^ej{b) + eij{b), (4) 

where 

/oo 
db\b'-b)j^^'[D{e)]ib'), (5) 
b 

eu{b) = (1 - ^){9'FF)T^'[^D^{k'm, (6) 

For r» = D or Di, 2?(fc2) = T4[V{z'^)], where Tn denotes a n-dimensional Fourier transform. 

With the above formalism, once one has inputs for the correlation functions D(z) and Di{z), the elementary 
amplitude in the momentum transfer space, Eq. (2), may, in principle, be evaluated through Eqs. (3 - 6). It is 
important to stress that, as constructed, this approach is intended for the high energy limit and small momentum 
transfer region, namely, s —> oo, where y/s is the cm. energy and q^ ^ C'(l) GeV^. 
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TABLE I: Central values of the fit parameters (without statistical errors) for the correlators (7) and (8) 0, 



parameters 


rUq.a — 0.01 


full QCD 


niq.a — 0.02 


quenched 
approximation 


An (fm^''') 


14.87 




31.04 


128.4 


ao 


0.71 




0.66 


0.69 


(fm-4) 


1.709 




4.102 


27.23 


tti 


0.45 




0.39 


0.46 


AA(fm) 


0.34 




0.29 


0.22 


\a (fm) 


4.4 




3.0 


0.43 



III. LATTICE PARAMETRIZATIONS 



The determination of the correlation functions through numerical simulation on a lattice is made by means of the 
cooling technique, a procedure that removes the effects of short-range fluctuations on large distance correlators. The 
numerical results with the associated error are usually parametrized with the functions [a, Q 

Diz) = Aoexp(-M) + ^exp(-M), (7) 

Aa \z\^ \a 



D,{z) = A, exp(-ii) + ^ exp(-U), (8) 

where A^i in the non-perturbative exponential terms is the correlation length of the gluon field strengths. Discussions 
on these choices, including the perturbative-like divergence at short distances, may be found in Refs. Q and 

These correlation functions were first determined in the quenched SU(3) theory and in the interval of physical 
distances (Euclidean space) between 0.1 and 1 fm Qlsj. After that, the effects of dynamical fermions have also been 
included (full QCD), for bare quark masses a.mq = 0.01 and a.niq = 0.02 (a is the lattice spacing) and physical 
distances between 0.3 and 0.9 fm 0. The above parametrizations are the same in all these cases and the numerical 
values of the parameters, obtained from references and 0, are displayed in Table 1. 



IV. RESULTS AND DISCUSSION 

With the procedure described in Sec. II (see Q for all the calculational details) , the elementary scattering amplitude 
in the momentum transfer space may be determined. The numerical results from full QCD, with a.mq = 0.01 and 
a.mq — 0.02, are shown in Fig. 1, together with our previous result in the quenched approximation. The normalized 
amplitudes are displayed in the region of high momentum transfer, up to 10 GeV^. 

We see that in all the cases the amplitudes decrease smoothly as the momentum transfer increases and they do not 
present any zeros in that region. The overall basic effect of the dynamical quarks is to originate a more rapid decrease 
of the amplitude, an effect that depends on the bare quark mass: the smaller the mass the faster the decrease. This 
behavior is associated with the correlation lengths and Aa, since they are the only parameters that decrease when 
going from full QCD (with increasing quark masses) to quenched approximation (see Table I). 

In order to obtain analytical expressions, suitable for investigating distinct contributions and also for phcnomcno- 
logical uses, we have parametrized these numerical points through a sum of exponentials in q^: 



m~h ■ 

By introducting a global uniform error of 0.5% in the numerical points, we fitted the data through the program 
CERN-Minuit. The results of the fits from full QCD with niq.a = 0.01 (approximation to the chiral limit) and in the 
quenched approximation are displayed in Table 2 and they are represented by the solid lines in Fig. 2, in the regions 
of large and small momentum transfer. The corresponding exponential components in each fit are shown in Fig. 3 
for the quenched approximation and in Fig. 4 for full QCD. 
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FIG. 1: Elementary amplitudes from full QCD and our previous result in quenched approximation. 



TABLE 11: Values of the fit parameters to the elementary amplitudes, Eq. (9), in the cases of quenched approximation and 
full QCD with m^.a = 0.01. 



Parameters: 


ai 

12 3 4 


A (GeV-2) 
12 3 4 


quenched 
full QCD 


0.03 - 0.69 0.28 
0.16 0.54 0.20 0.10 


55.0 - 0.57 0.09 
170.0 1.52 0.41 0.07 



An immediate conclusion from these results is the presence of an additional component in the case of full QCD. 
This is a central point that we shall discuss in certain detail in what follows. 

Let us start with the components with the highest slopes, which appear in both cases in the small momentum 
transfer region, namely, below ~ 0.1 GcV^ (Figs. 3 and 4). As mentioned before, the "real" lattice results 
correspond to sets of discrete theoretical points with errors, in a finite interval of physical distances, roughly 0.1 - 
1.0 fm. Therefore, the parametrizations (7) and (8) extrapolate this interval down and above. The highest physical 
distance reached in the simulations was 0.85 fm, which corresponds to ~ 0.24 GcV^. Therefore, we conclude that 
the components with the highest slopes in Figs 3 and 4 {q^ < 0.1 GeV-^) are associated with the extrapolations above 
the physical region with "real" lattice results. In the same manner, the components with the smallest slopes are 




2 4 6 8 10 0.2 0.4 0.6 0.8 1 

q' (GeV'^) q'' (GeV") 

FIG. 2: Fits to numerical points with paramctrization (9) in the cases of full QCD for rriq.a = 0.01 and quenched approximation, 
in the region of large (a) and small (b) momentum transfer. 



5 




q' (GeV^) (GeV'') 



FIG. 3: Exponential components of the fit at large (a) and small (b) momentum transfer in quenched approximation. 
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FIG. 4: Exponential components of the fit at large (a) and small (b) momentum transfer in full QCD for mq.a = 0.01. 



connected with extrapolations down the "real" lattice points in physical distances and they are the responsible for 
the amplitudes in high momentum region {q^ > Q — 7 GeV^ in Figs. 3 and 4); therefore, they are outsite the region 
of validity of the SVM, namely, < 0(1) GeV^. 

We conclude that in this context only the intermediate components, predominant in the interval, let us say, 0.5 < 
q"^ < 2.0 GeV^, can have physical meaning in the sense of being in agreement with the SVM and the "real" lattice 
results in both full QCD and quenched approximation. From Figs. 3 and 4 this "secure" region is characterized by 
only one component in the case of the quenched approximation and two components in full QCD. The transition 
between these two components occurs at ~ 1 GeV^ (Fig. 4 - left), a limit region for which the SVM is intended 
for. 



V. CONCLUSIONS AND FINAL REMARKS 

In this work we have obtained analytical paramctrizations for the quark-quark scattering amplitudes in a nonper- 
turbative QCD framework (SVM) and using as inputs the correlation functions, determined from numerical simulation 
on a lattice, in both quenched approximation and full QCD. The formalism is intended for small momentum transfer 
(g^ ^ 0{1) GeV^), asymptotic energies s ^ oo and physical distances between 0.1 and ~ 0.9 fm. As discussed 
in Sec. IV, these conditions put some restrictions in the physical interpretations and, therefore, in the practical 
phenomenological uses of these amplitudes. 

However, even under the above strictly conditions we can surely extract some novel results: (1) the amplitudes 



6 



decrease smoothly as the momentum transfer increases and they do not present zeros; (2) the decreasing is faster 
when going from quenched approximation to fuU QCD (with decreasing quark masses), and this effect is associated 
with the increase of the correlation lengths (A^ and A^); (3) the dynamical fermions generate two contributions in 
the region of small momentum transfer, which are of the same order at 1 GeV^ (only one contribution is present 
in the case of quenched approximation). 

We understand that result (3) may suggest some kind of change in the dynamics at the elementary level, near 
~ 1 GeV^ and at asymptotic energies. If that is true, some signal could be expected at the hadronic level. One 
possibility is that this effect can be associated with the position of the dip (or beginning of the "shoulder" ) in the 
hadronic (elastic) differential cross section data. The asymptotic condition embodied in our result indicates that ^ 
1 GeV^ seems in agreement with limit of the shrinkage of the diffraction peak, empirically verified when the energy 
increases in the region 23 GeV < ^/s < 1.8 TeV. 

At last it should be noted that if there is no new effect above the physical distances presently investigated on the 
lattice (^ 0.9 fm), the extrapolations can be considered as a good representation of the lattice results. In this case 
our analytical parametrizations may be usefuU inputs for phenomenological uses in the region of small momentum 
transfer and asymptotic energies. 
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